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Abstract 

The phenomenon of phase transitions in one-dimensional systems is discussed. 
Equihbrium systems are reviewed and some properties of an energy function which 
may allow phase transitions and phase ordering in one dimension are identified. We 
then give an overview of the one-dimensional phase transitions which have been studied 
in nonequilibrium systems. A particularly simple model, the zero-range process, for 
which the steady state is known exactly as a product measure, is discussed in some 
detail. Generalisations of the model, for which a product measure still holds, are also 
discussed. We analyse in detail a condensation phase transition in the model and 
show how conditions under which it may occur may be related to the existence of an 
effective long-range energy function. It is also shown that even when the conditions 
for condensation are not fulfilled one can still observe very sharp crossover behaviour 
and apparent condensation on a finite system. Although the zero-range process is not 
well known within the physics community, several nonequilibrium models have been 
proposed that are examples of a zero-range process, or closely related to it, and we 
review these applications here. 

1 Introduction 

In recent years the study of nonequilibrium systems has come to the fore in statistical 
mechanics. Basically, one can consider two types of nonequilibrium systems: those relaxing 
towards thermal equilibrium and those held far from thermal equilibrium e.g. by the system 
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being driven by some external field. In the present article we will be mainly concerned with 
the latter scenario. 

To be more specific we define our nonequilibrium systems as those evolving through a 
local stochastic dynamics which a priori does not obey detailed balance, at least not with 
respect to any 'reasonable' energy function. The question of what is a reasonable energy 
function is a moot point. One might propose that the energy contains only local interactions, 
or is extensive, or is written down according to some physical principles; but any answer 
to the question is subjective. However, the basic point is that the nonequilibrium system 
is defined by its dynamics without regard to any concept of energy and it is the dynamics 
which should seem reasonable or 'physical'. This is distinct from an equilibrium system 
where the energy function should be 'physical' and the dynamics is usually defined in an 
ad hoc way simply to guarantee that one obtains the Gibbs-Boltzmann weight with the 
specified energy. The easiest way to do this is to use the detailed balance condition. 

A natural way to construct a nonequilibrium steady state is to drive the system by forcing 
a current of some conserved quantity, for example energy or mass, through the system. Such 
systems are known as driven diffusive systems (DDS). The archetypal model was introduced 
Katz, Lebowitz and Spohn Basically it comprises a two dimensional Ising-like lattice gas 
evolving under conservative Kawasaki dynamics (spin exchange) and with a drive direction 
imposed. It has been shown that a continuous phase transition exists in the driven system, 
as is also the case in the undriven (Ising) system, but, most interestingly, one sees generic 
long range (power-law decay) correlations as opposed to the undriven systems where long- 
range correlations are only seen at criticality. Although exact results are not available for 
this systempl, it is often thought that generic power-law correlations are related to the 
existence of an effective long-range Hamiltonian for the system (see e.g. [0,3). 

More recently it has been realised that DDS in one dimension exhibit non-trivial be- 
haviour. The interest has been from a fundamental viewpoint but also in the context of 
applications such as interface growth ||^ and traffic flow modelling §]. Also it turns out that 
problems of transport with a single-flle constraint have long been of interest in biological 
contexts such as transport across membranes |^ and the kinetics of biopolymerisationP]. 

One intriguing feature of one-dimensional systems is the possibility of phase ordering 
and phase transitions. In recent years this possibility has begun to be explored and some 
examples are by now well studied. To appreciate the signiflcance one should recall the general 
dictum that in one-dimensional equilibrium systems phase ordering and phase transitions 
do not occur (except in the limit of zero-temperature, or with long range interaction — see 
section 2). In the one-dimensional nonequilibrium systems studied so far it appears that the 
presence of conserved quantities and an imposed drive are important in allowing ordering 
and phase transitions. However there still does not exist a general theoretical framework 
within which to understand the phenomena. 

The purpose of this article is twofold. Firstly I wish to give a broad overview of phase 
transitions and phase ordering in one dimension — this is carried out in Section In par- 
ticular, in Section ^?T| we discuss the conditions under which ordering and phase transitions 
may occur in equilibrium systems i. e. the requisite properties of the energy function to allow 
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such phenomena. Then in Section p.2| we catalogue some nonequihbrium, one dimensional 
systems which exhibit non-trivial phase behaviour 

The second purpose is to discuss a very simple class of microscopic models, the zero-range 
processes [|^, |lT|, which are presented in section For these models the steady state can 
be calculated exactly since it factorises into a product measure. There is some irony in the 
fact that the system has found widespread application in the modelling of nonequihbrium 
phenomena (see Section although the zero-range process was originally introduced by 
Spitzer|]10[ as a dynamics which could lead to Gibbs measures. In section 4 we discuss 
generalisations of the basic model which also have a product measure steady state. We 
show in Section |^ how the model can exhibit a phase transition, that we shall refer to as 
a condensation transition, which is analysed in some detail. We also discuss an interesting 
sharp crossover phenomenon whereby models, although not fulfilling the conditions for strict 
condensation and phase ordering, may often appear to be in a condensed phase on a finite 
system. The simplicity of the system allows us to explore the roles of a conserved quantity, 
the presence of a drive and effective long-range energy functions. Conclusions are drawn in 
Section ^. 



2 Phase transitions in one dimension 



2.1 One-dimensional equilibrium systems 

As mentioned above, it is received wisdom that in one-dimensional equilibrium systems 
phase transitions do not occur. In fact any careful statement of this requires a few caveats 
and, indeed, a general rigorous statement is hard to formulate (see for a discussion). 



Perhaps the best known argument is that of Landau and Lifshitz [|T^. For simplicity, 
consider a one dimensional lattice of L sites with two possible states, say A and B, for each 
site variable. Let us assume the ordered phases, where all sites take state A or all sites take 
state B, have the lowest energy, and assume a domain wall (a bond on the lattice which 
divides a region of A phase from that of B) costs a finite amount of energy e. Then n 
domain walls will cost energy ne but the entropic contribution to the free energy due to the 
number or ways of placing n walls on L sites ~ nT [ln(?7,/L) — 1] for 1 <C n -C L. Thus 
for any finite temperature a balance between energy and entropy ensures that the number 
of domain walls grows until it scales as L, that is, until the typical ordered domain size is 
finite. 

Note that this argument relies on a finite energy cost for domain walls, and short range 
interactions so that one may ignore the interaction energy of domain walls. Indeed, the Ising 
model with long-range interactions decaying with distance as J(r) ~ f'^^^" has been well 
studied [ff4[ and it has been demonstrated that the one- dimensional system orders at low 



temperatures for a < 1 . Also, of course we require non-zero temperature so that entropy 
comes into play, otherwise the two fully ordered states (ground states) would dominate the 
partition sum and the system would be frozen into them. 
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Another even simpler way of thinking of this is from a dynamical perspective. For 
a disordered state to order, domain walls must annihilate each other. However in one 
dimension no energy is gained by the two domain walls at opposite ends of a domain moving 
closer to one another; a domain always has two domain walls costing energy 2e no matter 
what its size is. Therefore there is no effective force to eliminate domains and the system is 
disordered. Again, this argument requires a short range interaction so that one can ignore 
the energy of interaction of domain walls above some finite distance. 

A more mathematical way of addressing the question of phase transition in Id is to use 
the transfer matrix technique For example, on a periodic one-dimensional homogeneous 



system of sites, the partition sum can be written as the trace of a product of A^ transfer 
matrices T : 

Z = Trace [T^] = (1) 

A 

where A are the eigenvalues of the transfer matrix. Now, since the transfer matrix is finite 
and the entries are all positive the Perron-Frobenius theorem tells us that the largest 
eigenvalue Amax is non-degenerate. Thus, there can be no crossing of the largest eigenvalue 
as we vary some control parameter. Consequently the free energy F oc \im.N^oo{\'^Z)/N = 
Amax is analytic and we have no phase transitions (which would be signalled by some non- 
analyticity of the free energy). 

Again, there are exceptions to this reasoning i.e. when the Perron-Frobenius theorem no 
longer applies. This can occur when the transfer matrix becomes infinite due, for example, 
either to long range interactions or when the local degree of freedom at each lattice site 
is not restricted to a finite number of states e.g. |T^. (An extreme instance of the latter 



case is when we are actually considering a two dimensional system!) Another case when the 
Perron-Frobenius theorem does not apply is when the transfer matrix becomes reducible i.e. 
when there exists components of that are zero for all values of A^. This can occur at 
zero temperature or when some interaction strengths are set to infinity, an example being 
the first order transition in the KDP model discussed in 



In this section we have discussed three arguments, presented here at different (low) 
levels of rigour, which all point to phase transitions in equilibrium one-dimensional systems 
only being possible in the case of long-range interactions, zero-temperature limit or infinite 
interaction energies, or unbounded local variable at a site. As we shall see the situation for 
nonequilibrium systems is less restrictive although some parallels can be drawn. 



2.2 One-dimensional nonequilibrium systems 

Here we give an overview of one-dimensional systems where phase transitions and phase 
ordering may occur. We focus our attention on hopping particle models that, despite their 
simplicity, offer a wide range of non-trivial behaviour. 

A simple one-dimensional model of a driven diffusive system is the asymmetric simple 
exclusion process (ASEP) . Here particles hop in a preferred direction on a one-dimensional 
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lattice with hard-core exclusion (at most one particle can be at any given site). Indicating 
the presence of a particle by a 1 and an empty site (hole) by the dynamics comprises the 
following exchanges at nearest neighbour sites 

10^01 with rate 1 

1 1 with rate q (2) 



The open system was studied by Krug||2^ and boundary induced phase transitions shown 
to be possible. Specifically one considers a lattice of N sites where at the left boundary site 
(site 1) a particle is introduced with rate a if that site is empty, and at the right boundary 
site (site A^) any particle present is removed with rate j3. Thus the dynamical processes at 
the boundaries are 

at site 1 — > 1 with rate a 

at site 1 ^ with rate (3 . (3) 

These boundary conditions force a steady state current of particles J through the system. 
Phase transitions occur when limTv^oo J exhibits non-analyticities. The steady state of this 



system was solved exactly for the totally asymmetric case pT| , p2| and more recently for 
the general q case p3| , p4| . When g < 1 the phase diagram comprises three phases: a high- 
density phase where the current is limited by a low exit rate (3 and takes the expression 
J = f3{l — q — j3) / {1 ~ q); a low-density phase where the current is limited by a low injection 
rate a and takes the expression J = a{l — q — a) / {1 — q)\ a maximal-current phase where 
both a,j3 > (1 — g)/2 and the current is J = (1 — g)/4. In the maximal current phase 
generic long-range correlations exist, an example being the decay of particle density from 
the left boundary to the bulk value 1/2 which is a power law ~ l/x^^^ where x is distance 
from the left boundary. 

Clearly the presence of a conserved quantity and a drive, leading to non-zero current J 
is crucial to the phase transition. Indeed, the qualitative phase diagram appears robust for 



stochastic one-dimensional driven systems ||2^. For the case of no bulk drive g = 1 p6|, 

or 'reverse bias' q > 1 the current vanishes with increasing system size and there are no 

boundary-induced phase transitions. 

The model has been generalised to two oppositely moving species of particle: one species 
is injected at the left, moves rightwards and exits at the right; the other species is injected 
at the right, moves leftwards and exits at the left [Q. Spontaneous symmetry breaking 
has been shown to occur, whereby for low exit rates {j3) the lattice is dominated by one 
of the species at any given time. In the low (3 limit the mean flip time between the two 
symmetry-related states has been calculated analytically and shown to diverge exponentially 



with system size 29 



In these models the open boundaries can be thought of as inhomogeneities where the 
order parameter (particle density) is not conserved. Inhomogeneities which conserve the 
order parameter can be considered on a periodic system. Indeed a single defect bond on 
the lattice (through which particles hop more slowly) is sufficient to cause the system to 
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separate into two macroscopic regions of different densities : a high density region which 



can be thought of as a traffic jam behind the defect and a low density region in front of the 
defect. Here the presence of the drive appears necessary for the defect to induce the phase 
separation. 

Moving defects {i.e. particles with dynamics different from that of the others) have also 
been considered and exact solutions obtained |31|, ^ ^ In the model studied 
33| , 0, varying the rate at which the defect particle hops forward, denoted a, and 
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the rate at which it is overtaken and exchanges places with normal particles, denoted /3, 
produces a phase diagram closely related to the open boundary problem. Moreover for low 
f3 and high a there is a phase where the defect particle induces phase separation between a 
high density region behind it and a low density region in front of it. 

For some of the examples discussed so far the steady state has been solved exactly by 
constructing a matrix product which is reviewed in This reveals that the steady state 



weights are very complicated functions of the particle number and positions. It does not 
appear easy to relate this to any concept of an energy function. Indeed, it has been shown 
that a matrix product state is non-Gibbsian . 



A natural question to ask is whether systems related to the hopping particle models 
described so far, but without inhomogeneities, can exhibit phase ordering. A very simple 
model was introduced in PB| comprising three species of conserved particles, amongst which 



all possible exchanges are allowed. However a key feature is that the dynamics has a cyclic 
symmetry. To be specific let each site of a one-dimensional periodic lattice be occupied 
either by an A,B ot C particle (there are no holes in this model). The dynamical exchanges 
are 



AB - 


^ B A 


with rate 


q 


B A - 


^ AB 


with rate 


1 


B C - 


CB 


with rate 


q 


C B - 


B C 


with rate 


1 


CA - 


^ AC 


with rate 




AC - 


CA 


with rate 


1 



(4) 

and we will take g < 1. For example, the hopping of an A particle is biased to the right 
when it is an environment of Cs and it is biased to the left when it is in an environment of 
Bs. 

The phase separation observed in the model is rather easy to understand: if the system 
has separated into a domain of As, followed by a domain of 5s, followed by a domain of Cs 
(in that order), then the domain walls that are present AB., BC, CA are all stable objects. 
This is clear from (^) since, for example, any A particles which penetrate the B domain 
will be driven backwards by the dynamics. On the other hand BA, CB or AC walls are all 
unstable objects and would be quickly eliminated by the dynamics. 

In the special case of exactly equal numbers of A, B and C particles it was shown that 
the model actually obeys detailed balance with respect to a long range asymmetric, energy 
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function. In fact the energy is non-extensive in the sense that most configurations have 
energies of order A^^ . The partition sum was calculated in the large limit (with q fixed) 
and shown to depend linearly on N . This reflects the fact that the phase separation is into 
three pure domains and the partition sum is dominated by the equivalent translations 
of the structure comprising three pure domains. When the numbers of particles are not 
identical, detailed balance does not hold but the phase separation into pure domains remains. 
Similar behaviour has been found in other related models with conserving dynamics ^ 



Another interesting model is where phase separation occurs on a quasi-one-dimensional 
system [2 x N sites) but not on a strictly one-dimensional system |iT]. It should also be 



mentioned that systems with a cyclic symmetry but with non-conserving dynamics have 



been studied and shown to order into a frozen state |H2 



Any discussion of nonequilibrium phase transitions is not complete without mentioning 
the most well known class, that of directed percolation. Various models are reviewed else- 
where in this volume [^] so here I just sketch the basic behaviour by referring to a particular 



model, the contact process p4| , Each site of a lattice is either empty or contains a par- 
ticle. Particles are annihilated with rate 1 and particles are created at empty sites with rate 
nX/2 where n is the number of occupied nearest neighbours of the site [n = 0, 1,2). Note 
that the 'inactive state' where all sites are empty is an absorbing state. Above a critical 
value of A there is a finite probability that starting from a single particle on an infinite 
lattice, the system will remain active as t — >■ oo. This phase transition has well-studied as- 
sociated critical exponents and scaling behaviour. Moreover it appears to be a universality 
class in the sense that the same exponents are found in all systems, with the same symmetry 
and conservation laws, exhibiting a phase transition from an absorbing inactive state to an 



active state 45 



However as described so far the contact process is distinct from the other hopping particle 
models discussed in that on any finite lattice the absorbing state is reached in a finite time 
and is therefore the steady state. The active state only becomes available as a steady state on 
an infinite system. We mention briefly that it is in fact possible to define hopping particle 
models, similar in spirit to the nonequilibrium models discussed in previous paragraphs, 
that exhibit phase transitions connected with directed percolation. These models can have 
non-conserved order parameter |^B| or conserved order parameter Although there are 



no absorbing states in these models, they have the common feature of certain microscopic 
processes being forbidden. 

A final class of transitions in one-dimensional hopping particle models is that involving 
spatial condensation, whereby a finite fraction of the particles condenses onto the same 
site. Examples include the appearance of a large aggregate in models of aggregation and 
fragmentation p8[ and the emergence of a single flock in dynamical models of flocking |^9|, ^ 



In the Section 5 we shall examine a very simple example of a condensation transition which 
occurs in the zero-range process and see how it is related to a defect induced transition. 
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3 The zero-range process 



The zero-range process was first introduced into the mathematical hterature as an example 
of interacting Markov processes [|l^]. Since then the mathematical achievements have been 
to prove existence theorems, invariant measures and hydrodynamic limits [0, 

It is not widely appreciated that the zero-range process has many physical applications; 
moreover it has often appeared incognito in a wide range of different contexts. Examples 
include the repton model of polymer dynamics with periodic boundary conditions [Q; a 
model of sandpile dynamics [Q; the backgammon model |^4| for glassy dynamics due to 
entropic barriers; the drop-push model for the dynamics of a fluid moving through backbends 
in a porous medium [^; microscopic models of step flow growth |^6|, |5^ and a bosonic lattice 
We shall discuss some of these in the sequel. The zero-range process is also closely 

m 



58 



gas 

related to the more widely known asymmetric exclusion process 
below. 
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as we shall describe 



3.1 Model definition 

In general one can consider the zero-range process on a lattice of arbitrary dimension, and of 
(countably) infinite or finite number of sites. Initially Spitzer [|T^] considered a finite number 
of sites. However, subsequently most mathematical works tackle the invariant measure on 
an infinite system |TT]. For our purposes, it is most convenient to consider a finite system. 



compute the steady state and only then take the limit of an infinite system. Note that the 
steady state t — oo and the infinite volume limit do not necessarily commute e.g. on an 
infinite system the invariant measure (steady state) is not necessarily unique. 

We consider a one- dimensional finite lattice of M sites with sites labelled n = 1 . . . M 
and periodic boundary conditions. Each site can hold an integer number of indistinguishable 
particles. The configuration of the system is specified by the occupation numbers of each 
site /i. The total number of particles is denoted by L and is conserved under the dynamics. 
The dynamics of the system is given by the rates at which a particle leaves a site fi (one 
can think of it as the topmost particle — see Figure la). As our first example we assume 
it moves to the left nearest neighbour site /i— 1. The hopping rates u{n) are a function of 
n the number of particles at the site of departure. Some particular cases are: if u{n) = n 
then the dynamics of each particle is independent of the others; if u{n) = const for n > 
then the rate at which a particle leaves a site is unaffected by the number of particles at the 
site (as long as it is greater than zero). It is helpful to think of performing a Monte-Carlo 
simulation: in the u{n) = n case at each update a particle would be picked at random and 
moved to its nearest neighbour site; in the u{n) = constant case a site would be picked at 
random and a single particle moved to the nearest neighbour site. 



A possible source of confusion in the definition of the model is that in []TU[ and some 
other papers the hop rates u{n) are defined as the hop rate j)er particle at a site; thus u{n) 
in those works are 1/n of the u{n) defined here. 
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The important attribute of the zero-range process is that it yields a steady state described 
by a product measure. By this it is meant that the steady state probabihty P({n^}) of 
finding the system in configuration {^1,77,2 . . - Um} is given by a product of factors f{n^) 
often referred to as marginals 

1 ^ 

= 1(^7)11 

Here the normalisation Z{M, L) is introduced so that the sum of the probabilities for all 
configurations, with the correct number of particles L, is one. We shall explore later in 
Section ^ the interesting possibilities afforded by the form (|^). 

In the basic model described above, f{n) is given by 

n ^ 

u{m) 

m=l ^ ' 

= 1 for n = (6) 



Note that f{n) is defined only up to a multiplicative constant and we could have included 
a factor in (H). Later this factor reappears as a fugacity in Section ^. 

The proof of (^,|^) is, happily, straightforward. One simply considers the stationarity 
condition on the probability of a configuration (probability current out of the configuration 
due to hops is equal to probability current into the configuration due to hops): 

9{n^)u{n^)P{ni . . . . . . n^) = ^ e{n^)u{n^+i + l)P{ni . . . ...ul). 

(7) 

We have included the Heaviside function to highlight that it is the sites with n > 1 that 
allow exit from the configuration (Ihs of (0)) but also allow entry to the configuration (rhs 
of (1^). Equating the terms fj, with fj, > 1 and cancelling common factors assuming (|^), 
results (for n^> 1) in 

u{n^)f{n^_i)f{n^) = M(n^+i + l)/(n^ - l)/(^^^+i + 1) (8) 

This equality can be recast as 

u{n^) ,/^^^\. = u{n^+i + + = constant (9) 



Setting the constant equal to unity implies 



/K - 1) 



fM = ^ / (10) 



and iterating ( |10|) leads to (||) where we have chosen /(O) = 1. 
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12 3 4 

particle 

Figure 1: Equivalence of zero range process and asymmetric exclusion process. 

3.2 Relation to the asymmetric exclusion process 

There exists an exact mapping from a zero-range process to an asymmetric exclusion process. 
This is illustrated in Figure 1. The idea is to consider the particles of the zero-range 
process as the holes (empty sites) of the exclusion process. Then the sites of the zero-range 
process become the moving particles of the exclusion process. This is possible because of 
the preservation of the order of particles under the simple exclusion dynamics. Note that 
in the exclusion process we have M particles hopping on a lattice of M + L sites 

An interesting feature of the mapping is that it converts a model where the local degree of 
freedom can take unbounded values (particle number in the zero-range process) to a model 
where the local site variable is restricted to two values. On the other hand, a hopping rate 
u{m) which is dependent on m corresponds to a hopping rate in the exclusion process which 
depends on the gap to the particle in front. So in principle the particles can feel each other's 
presence and it is possible to have a long-range interaction. 

4 Generalisations 

We now show how the totally asymmetric, homogeneous zero-range process we have consid- 
ered so far may be generalised yet retain steady states of a similar form to d^J^). 
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4.1 Inhomogeneous system 



First we consider an inhomogeneous system by which we mean the hopping rates are site 
dependent: the hopping rate out of site ^ when it contains particles is u^{n^). It is easy 
to check that the steady state is simply modified to 

1 ^ 

p^M) = ziMjr)IlfM- (11) 

where are given by 

n ^ 

/m(^) = n — ^ - 1 

J- J- 1/ I m 1 



m=l ' 



= 1 for n = (12) 

The proof is identical to that given above for the homogeneous case, with the trivial re- 
placement of M(n^) by u^{n^) 



4.2 Discrete Time Dynamics 



A further generalisation is to the case of discrete time dynamics, 
in 



This has been studied 



59| in the context of a disordered asymmetric exclusion process. Here we translate the 



results into the zero-range process. Rather than processes occurring with a rate, time is 
counted in discrete steps and at each time step events occur with certain probabilities. 



In the case of Parallel Dynamics, at each time-step all sites are updated. One particle 
from each site yU is moved to the left, each with probability p^{n^) where is the number 
of particles at the site before the update. Note that the particles move simultaneously and 
particles do not move more than one site. 



It turns out that the steady state again has the form ([TT]) . It was shown in that 
ff,{n) = l-p^(l) for n = 

1 -P/x(l) TT 1 -Pt^i^^ 



1-M 



n 



n 

m=l 



for n > 



m] 



(13) 



To recover the continuous time dynamics we can call the interval between time steps dt 
and let p^{n^) = u^{n^)dt. Then continuous time dynamics is given by the limit dt ^ Q 
and, to within a constant factor dt", (|13D recovers ([T^ . In this way one can interpolate 
between discrete time, parallel dynamics and continuous time dynamics. 

In ordered sequential updating schemes were also considered. These are discrete 
time updating schemes were one site is updated at a time, but the sites are updated in a 
fixed order. The steady states for the forwards and backwards updating sequences were 
derived and it turns out they too have the form (|r^) with taking an expression related 
to the parallel case (|T^). 
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4.3 Arbitrary Network 



In the original paper paper of Spitzer|TO| some more general versions of the zero range 
process were considered. Here we discuss one interesting case which serves to generalise the 
(totally asymmetric) zero range process defined above to a process on a more general lattice 
or for any finite collection of points with a prescribed transition matrix for the dynamics of 
a single particle 



In this case the rate of hopping of a particle at site fi containing particles is equal 
to u^{n^) and the probability that a particle leaving site /i will move to site z/ is denoted 
Wlfi-^u). Thus the probability that in time dt a particle at fi moves to u is 

u^in^)Wifx^u) dt . (14) 

Note that the probabilities Wlfi^u) define a stochastic matrix for a single particle moving 
on a finite collection of M sites and we take 

J2wifi^u) = l. (15) 

V 

so that probability is conserved. We refer to the collection of points together with the 
prescribed transition matrix Wi^^-^v) as the network. 

We assume that the transition matrix is irreducible [i.e. the particle can pass from any 
given point to any other after sufficient time and the system is ergodic) so that we have a 
unique steady state probability for the single particle problem: 

= ^s^W^(/i^z/) . (16) 



We now show that the steady state for the many-particle problem defined above is given 
by ( pAD where now f^{n) is given by 

n 

uin) = n 

m=l 

= 1 for n = (17) 

The proof is again a straightforward generalisation of that of Section |3T| . Equation (|^) is 
modified to 

^ 9{n^)u^{n^)P{ni . . . n^) = ^ ^ e{n^)p{u-*n)u^{n^+l)P{ni . . . n„+l . . . ...Ul). 

(18) 

Equating the terms /x on each side of (|18]), assuming ([TT| ) and cancelling common factors 
yields 



uJm) 



for n > 1 
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Inserting ( |I7D leads to the condition 



which is the same as the single particle steady state condition (0). 



(20) 



A simple case considered by Spitzer is when Wlfi-^u) is a doubly stochastic matrix 
which is defined by the property 



(21) 



Equations (^T]) and (jl^) then imply that the single particle problem has a homogeneous 
steady state = constant. 

Let us also discuss an example where the single particle problem has an inhomogeneous 
steady state. We consider a one- dimensional lattice where hops to the left and right neigh- 
bours are allowed but with probabilities that depend on the site. Thus, we may write 



W{ii—>-i') = for z/ = /i + 1 

= 1 — for z/ = /i — 1 
= otherwise . 



(22) 
(23) 
(24) 



The steady state of the single particle problem (random walker on a disordered one dimen- 
sional lattice 16011) 



= (1 - g^+i)s^+i + g^_is^_i 



can be solved and one obtains 



'M-l i 



(25) 



(26) 



This network is relevant in the disordered one-dimensional exclusion process studied 



m 



61, 62, 53]. The sites in the present model correspond to the particles in the exclusion 



process which each have their own forward and backward hopping rates. Another, particular 
instance of this network occurs in [^], where a repton model of gel electrophoresis is studied 
in the case of periodic boundary conditions (see Section 

In special cases of the zero-range process detailed balance may hold. The condition for 
this is 

u^{n^)W{iJi-^u)P{ni . . .til) = Uu{n^ + l)W{u-^fi)P{ni . . . n^+l . . . n^-l . . .ul) ■ (27) 
Substituting the form ([TT|JT7| ) leads to the condition 

s^W{^l^v) = s,W{v^^^) (28) 



13 



which is just the detailed balance condition for the single particle problem. 

An interesting consequence of the form of the steady state ([T7|) is that it allows one 
to relate an arbitrary zero-range process to a model obeying detailed balance. The idea 
is that if detailed balance doesn't hold, we can always define a new zero-range process (to 
be denoted by a prime) with the same steady state, but with a different dynamics obeying 
detailed balance. To do this, we solve the single particle problem (|T6|) for the original 
model to obtain s^. For any collection of points we can always define a new single particle 
transition matrix W'lfi^u) that satisfies detailed balance with respect to a homogeneous 
steady state (s^ = constant). The new model is defined by a new set of hopping rates 
u'^{m) = Ufj^{m)s'^/Sfj_ together with the new transition matrix W'lfi^u). It is easy to check 
from ([T7|) that the new model has the same steady state as the original. 

Thus, within the realm of zero-range processes, to the steady state of any nonequilibrium 
model we can always identify a model satisfying detailed balance and therefore an energy 
function. Of course, although the steady states are the same, there is no reason for the 
dynamical behaviour of the two systems to be related. To clarify this point we will discuss 
a simple example in Section 

The marginals (|1^) have the interesting structure of being a product of a term (s^)" 
that depends on the nature of the network and a term involving the product of Uy{m) which 
reflects the interactions at the site. The network can represent an arbitrary dimensional 
lattice or the effects of disorder, the only difficulty to surmount in obtaining the steady 
state is the solution of the single particle problem. 



5 Condensation Transitions 

We now proceed to analyse the steady states of form ([TT|) and the condensation transition 
that may occur. The important quantity to consider is the normalisation Z{M, L) as it 
plays the role of the partition sum. The normalisation is defined through the condition 

M 

Z{M,L)= '^(E^^-^)n/MK) (29) 

where the 6 function enforces the constraint of L particles. The normalisation may be 
considered as the analogue of a canonical partition function of a thermodynamic system. 

We define the 'speed' v as the average hopping rate out of a site 

M 



^ Z(M L) ^ ^(J2^>^-LHn,)f[f^{n^) 

' ni,n2...nM M M=l 

_ Z{M,L-1) 

~ Z{M,L) ^"^"^ 

where we have used (|TT||T^). Note that (^) tells us that the speed is independent of site and 
thus may be considered a conserved quantity in the steady state of the system. In the totally 
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asymmetric system considered in Section the speed is equal to the current of particles 
flowing between neighbouring sites and is clearly a conserved quantity in the steady state. 
More generally, however, the speed is not equal to the current and the fact that the speed 
is a conserved quantity is not a priori obvious. For example, in a system obeying detailed 
balance the (net) current is zero, but the speed as defined above remains finite. The speed 



is a ratio of partition functions of different system sizes (|30D and corresponds to a fugacity, 
as we shall see below. 

We will consider also the probability distribution of the number of particles at a site, 
taken here to be site 1 

1 ^ 
Pi(ni=x) = ^ ^(^ + ^2 + ■■■nM-L)fi{x)Ylf^{n^) 

' ' n2...nM M=2 

^ Z{M-l,L-x) 
= ^^("^ ZiM,L) 

(where Z{M—1, L—x) is the partition function for a system with site 1 removed). In general 
the probability distribution is site dependent but for a homogeneous system (/^ independent 
of /i) it will be the same for all sites. 

We now use the integral representation of the delta function to write the partition func- 
tion as 

M 



(32) 



where 

oo 

F,{z) = 5^ /,(m) . (33) 



m=Q 



For large M, L (|32|) is dominated by the saddle point of the integral and the value of z at 



the saddle point is the fugacity. The equation for the saddle point reduces to 

± = ±y^l^FJz) (34) 
which, defining = L/M, can be written as 

In the thermodynamic limit, 

M ^ oo with L = (j)M , (36) 

where the density (p is held fixed, the question is whether a valid saddle point value of z 
can be found from (|35|). We expect that for low (p the saddle point is valid but, as we shall 
discuss, there exists a maximum value of z and if at this maximum value the rhs of (^5]) is 
finite, then for large (^) cannot be satisfied. We now consider separately, and in more 
detail, how condensation may occur in the inhomogeneous and the homogeneous case. 
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5.1 Inhomogeneous case 



In general, the inhomogeneous case i.e. where F^{z) depends on the site yU through (|T3), is 
difficult to analyse. Here we would just like to give an idea of how a condensation transition 
may occur by discussing a simple example. We then go on to analyse perhaps the simplest 
example of a condensation transition: a single inhomogeneous site [Q. 

First we take the general model discussed in Section 4.3 and set u^{m) = u^foi m > i. e. 
the hopping rate does not depend on the number of particles at a site. We consider doubly 
stochastic transition matrices Wlfi-^u) (see Eq. ^) so that we may take s(/i) = constant 
and without loss of generality we set the constant equal to one. For the moment we do not 
specify further the transition matrix; later we will discuss two specific examples one obeying 
detailed balance and one not. Under these conditions is given by 

Un) - (37) 

and the probability of occupancies {ni, n2, . . . , n^} is 

1 *^ / 1 X*^" 

nM})=^(^n(-) . (38) 

The mapping to an ideal Bose gas is evident: the L particles of the zero-range process are 
viewed as Bosons which may reside in M states with energies E^^ determined by the site 
hopping rates: exp{—PE^) = Thus the ground state corresponds to the site with the 

lowest hopping rate. The normalisation Z{M, L) is equivalent to the canonical partition 
function of the Bose gas. We can sum the geometric series (|33|) to obtain and F'^ then 
taking the large M limit allows the sum over /i to be written as an integral 

duViu) ^— (39) 
u — z 

n 

where V[u) is the probability distribution of site hopping rates with Umin the lowest possible 
site hopping rate. Interpreting V{u) as a density of states, equation (|39D corresponds to the 
condition that in the grand canonical ensemble of an ideal Bose gas the number Bosons per 



state is 0. The theory of Bose condensation |64| tells us that when certain conditions on the 



density of low energy states pertain we can have a condensation transition. Then (pSD can 
no longer be satisfied and we have a condensation of particles into the ground state, which 
is here the site with the slowest hopping rate. This case is discussed further, in the context 
of an asymmetric exclusion process on an infinite system, by J. Krug in this volume ||65|| . 



We now turn to the simplest case of an inhomogeneous system: site 1 has Ui = p while 
the other M — 1 particles have hopping rates = 1 when /i > 1. It is easy to see that (0) 
simplifies to 
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In this case the normahsation Z{M, L) is easy to calculate combinatorially: 



Z{M,L) = Yl ^(E^A^-^)^"'^' 

p-^^ , (41) 



n\,n2...nM M 
L 



L+M-ni-2 
M-2 



= E 

ni=0 

yielding an exact expression for the speed through (|30|). In the thermodynamic limit the 
sum ( PI) is dominated by rii ~ 0{1) for (p < p/{l —p) and rii ~ C^(-^^) for (j) > p/{l —p) 
and it can be shown that 

for '^<7^ ^l/lTi^;;3jT3^ and. ^ 1-p (42) 



l-p V ' ^ I M ] i+4> p-(P(i-p) 

for </) > Z{M, L) ^ p-^ (1 - j9)-(^^-i) and v (f)/{l + 0) (43) 

1 — p 

In the high density phase, defined by (^) we have a condensate since the average number 
of particles at site 1 is {ni)/L = — p/(l — p). In the low density phase (^) the particles 
are evenly spread between all sites and we will refer to it as the homogeneous phase. 

We now discuss two models which both have this steady state: a driven system and 
a system obeying detailed balance. This provides an illustration of the idea discussed in 
Section 4.3 whereby a zero range process not obeying detailed can be related to one obeying 
detailed balance. 

First we take the totally asymmetric model so that particles move to the site to the left: 
the transition matrix is 



So this model is similar to that discussed in Section |3.1| , and a mapping to a totally asymmet- 
ric exclusion process can be made in the same way as Section |3.2| . The equivalent exclusion 
process is illustrated in Figure 2. We see that the equivalent exclusion process is system 
of hard-core particles hopping to the right, one particle being slower than the rest. The 
interpretation of the two phases within the context of the exclusion process is that in the 
condensed phase (for the exclusion process a low density of particles) a 'traffic jam' forms 
behind the slow particle and the slow particle has a finite fraction of the lattice as 'empty 
road' ahead. Whereas in the homogeneous phase (a high density of particle for the exclusion 
process) the particles are roughly evenly spaced. 

On the other hand one may consider the case where the one particle problem is a sym- 
metric random walk so that the system obeys detailed balance. The transition matrix is 
given by 

When we map this system to a simple exclusion process we see from Figure 2 that we have 
a system of particles, the bulk of which perform a symmetric exclusion dynamics but with 
two adjacent asymmetric particles: the left one biased to the left and the right one biased 
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1/2 1/2 1/2 p/2 p/2 1/2 1/2 1/2 



Figure 2: A totally asymmetric exclusion model (upper) and a model with zero current 
(lower) that have equivalent steady states leading to (^if) . 



to the right. In the condensed phase the gap between these particles diverges. Previously a 
single asymmetric particle in a sea of symmetric particles has been studied but in that 
case there is no phase transition. At first it seems that we have found a counterexample to 
the received wisdom that no phase transition should occur in an equilibrium system, since 
we have a condensation transition in a model with local dynamics obeying detailed balance. 
Inferring an energy function from the steady state ( P0| ) through the following equation 

exp [—{PE)] = exp [— (x2 — xi) \np] 

reveals that our effective energy increases linearly with distance X2 — xi between the two 
asymmetric particles. Therefore the energy is 'unphysical' in that it has very long range in- 
teractions. Thus the phase transition can be rationalised within the categories of exceptions 
discussed in Section 

We have seen that this simplest example of a condensation transition (a single inhomo- 
geneous site in the zero range process) is exhibited both in a driven model and also in a 
model obeying detailed balance but with long-range energy function. Again it should be 
stressed that although the steady states of these two models are equivalent, the dynamical 
properties should be very different. For example in the homogeneous phase of the driven 
model we expect asymmetric exclusion like behaviour and the dynamic exponent should 
be 3/2 implying relaxation times of M^^^ on a finite system. However in the homogeneous 
phase of the model obeying detailed balance we expect symmetric exclusion like behaviour 



and the djTiamic exponent to be 2 implying relaxation times of M [37 



5.2 Homogeneous case 



We now consider the homogeneous zero-range process where in (14) the hopping rates are 



site independent and the single particle problem (^) has a homogeneous steady state = 1 
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|68| . A similar analysis has been carried out in the context of balls-in-boxes and branched 
polymer models pP| . 

In the present case, (^) is independent of /i and reads 



oo n 



n=0 m=l 



u[m] 



(44) 



The fugacity z must be chosen so that F converges or else we could not have performed 
(P3|). Therefore z is restricted io z < (3 where we define (3 to be the radius of convergence of 
F{z). From (|4^ ) we see that (3 is the limiting value value of the u{m) i.e. the limiting value 
of the hopping rate out of a site for a large number of particles at a site. We interpret (^) 
as giving a relation between the density of holes (number of holes per site) and the fugacity 
z. The saddle point condition ( ^5]) becomes 

= -fW ^ ^ 

Given that the rhs of (^Sj) is a monotonically increasing function of z (which is not difficult 
to prove) we deduce that density of particle increases with fugacity. However if at z = /?, the 
maximum allowed value of z, the rhs of (^) is still finite then one can no longer solve for the 
density and one must have a condensation transition. Physically, the condensation would 
correspond to a spontaneous symmetry breaking where one of the sites is spontaneously 
selected to hold a finite fraction of the particles. 



Thus, for condensation to occur {i.e. when is large enough for (^5|) not to have a 
solution for the allowed values of z) we require 

F'(z) 

lim —fl < oo . (46) 

z^l3 F{z) ^ ' 

We now assume that u{n) decreases uniformly to (3 in the large n limit as 

u{n) = (3{l + an)) (47) 
where C(^) is a monotonically decreasing function. Analysis of the series 

oo ( n 

F{(3) = J]expi-5^1n[l + C(m)] 

n=0 L m=l 

oo n 

F'{(3) = J]nexpi-5^1n[l + C(m)]}> (48) 



m=l 



reveals that the condition for condensation is simply that F'{(3) is finite and this occurs if 
u{n) decays to P more slowly than /3(1 + 2/?t,). (This is easiest to see by expanding In [1 + C] 
and approximating the sum over m by an integral in 



In order to fit this result into the picture of section Section |2.1| one can argue that 



since the condensate has an extensive number of particles at a site, the local site variable is 
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unbounded. Therefore the 'no phase transition rule' does not apply. One also gains insight 
by translating the results into the language of the simple exclusion process. In this context 
we can have condensation if the hop rate of a particle into a gap of size n decays as /3(l+2/n) 
therefore there is an effective long range interaction. 



5.3 Sharp crossover phenomena 

Having discussed the case where a true phase transition occurs we now consider a homoge- 
neous example where, although there is no strict condensation transition, some interesting 



crossover phenomena occur ||68 
Consider 



u{n) = 1 for n < r 
u{n) = (3 for n > r . 



(49) 
(50) 



One can interpret these hop rates as meaning that a site only distinguishes whether it 
contains greater than r particles. When we use the mapping of section 3.2 to a totally 
asymmetric exclusion process r becomes the range of the interaction in the sense that it is 
the number of sites ahead upto which a particle in the exclusion process distinguishes. 



When these hopping rates are inserted in (P^ one obtains 

r— 1 CO 

F(^)=5^^" + 5^^«r 



(51) 



n=0 



Performing the geometric series readily yields 



F{z) 
zF\z) 



+ 



I- z^ 



rz 



r-l 



+ 



rz 



r-l 



{l-zf \-z {(5-zf {(5-z) 



(52) 
(53) 



Then we find the condition (|45|) can be written after a little algebra as 



(/5-;^)2[0-z(l + 



z^\\ - P) [{l+P-2z)z - ^[l-z)iP-z)] 
+rz'Xl-z){(3-z){l-(3) 



Therefore for large r we find the solutions 



for 
for 
for 



< 



> 



1-/? 

£_ 
1-/3 



1-/3 



1+0 ■ \1+<I>J (!+</>) (<t>-0(,l+4>)) 

z^[3- /3(i+'-)/3(l - /?) 



(54) 

(55) 
(56) 
(57) 
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Thus we see as r—>-oo we have two phases: a high density phase (j) > /?/(! — /?) where 
the speed is (3 and a low density phase where the speed is (/)/(! + 0). In fact these phases 
correspond exactly to those of the single defect problem discussed in the previous subsection 
^) with j3 playing the role of p. For finite r, z is actually a smoothly increasing function 
of but we see from ( p5| , p7D that the curve sharpens as r increases. This is illustrated in 
Figure 3 where the numerical solution to (|^) is plotted is plotted for r = 10, 20, 30. One 
sees a dramatic sharpening as r increases leading to a sharp crossover between a low density 
and high density regime. 



0.5 



Z 



0.45 



0.4 








Figure 3: Solutions to (|5^) for (3 = 0.5 and r = 10, 20, 30 (increasing in sharpness of curve). 



In order to see the effects of this sharp crossover it is interesting to consider the particle 
number probability distribution ( pi]) which for this system is site independent and given by 



P{x) ~ for X < r 

~ {z/py for X >r 



(58) 
(59) 



One can think of this as a sum of two distributions, one for poorly occupied sites x < r and 
one for well occupied sites x > r. When > /?/(l — /?) the probability distribution for large 

X > r goes as 



P(x) 



exp 



)/2. 



f3 



(</.-/?(l + 0))V2 



(60) 
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so that the typical occupancy of well-occupied sites goes as /?~^/^. Taking, for example 
(3 = 0.1 and r = 10 leads to a typical occupancy of ~ 10^. Therefore to simulate the model 
one requires a number of particles very much larger than this! If care is not taken to do this, 
and the total number of particles in the system is comparable to the typical occupancy, one 
would have an apparent condensate on a finite system. 

An example of this phenomenon was studied recently within the context of a 'bus route 
model' [0 . There the underlying motivation was to consider how a non-conserved quantity 
could mediate an effective long-range interaction amongst a conserved quantity in a driven 
system with a strictly local dynamical rule. The model considered was defined on a Id 
lattice. Each site (bus-stop) is either empty, contains a bus (a conserved particle) or contains 
a passenger (non-conserved quantity). The dynamical processes are that passengers arrives 
at an empty site with rate A; a bus moves forward to the next stop with rate 1 if that stop is 
empty; if the next stop contains passengers the bus moves forward with rate f3 and removes 
the passengers. Since the buses are conserved, there is a well defined steady state average 
speed V. This fact can be used to integrate out the non-conserved quantity (passengers) 
within a mean-field approximation. The idea is that a bus stop, next to bus 1 say, will last 
have been visited by a bus (bus 2) a mean time ago oin/v where n = X2 — a;i is the distance 
between bus 2 and bus 1. Therefore the mean-field probability that the site next to bus 1 
is not occupied by a passenger is exp(— An/f). From this probability an effective hopping 
rate for a bus into a gap of size n is obtained by averaging the two possible hop rates 1,(3: 

M(n) =/3 + (l-/5)exp(-An/f) . (61) 

We can now see that this mean-field approximation to the bus-route model is equivalent 
to a homogeneous zero-range process as discussed earlier in this section. Since u{n) decays 
exponentially, with decay length r = v/X, the condition for a strict phase transition is not 
met. It is reasonable to believe that the system behaves in a similar way to the system with 
a finite 'range' r discussed in Section p.3| . Since r can be made arbitrarily large as A— i>0, 
on any finite system an apparent condensation will be seen. In the bus route problem this 
corresponds to the universally irritating situation of all the buses on the route arriving at 
once. 



6 Some further applications 

As mentioned earlier the zero-range process and related models have appeared several times 
in the modelling of nonequilibrium phenomena. Here we briefly discuss a few of these 
instances to illustrate the ubiquity of the basic model. 



In [|5^ models of sandpile dynamics are considered. A zero range process is used to model 
the toppling of sand on a one-dimensional lattice; specifically the system is homogeneous and 
the occupation number of a site becomes the height of sand (h) at that site. The hopping 
rates are set as u{l) = 1 and u{h) = A for h > 1, with the transition matrix a symmetric 
random walk, and the limit of large A considered. This limit means that a particle (grain of 
sand) keeps moving until it finds an unoccupied site, thus a hopping event may play the role 
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of an avalanche. (Although in terms of sandpiles and self-organised criticality this model is 
rather trivial, it did serve to investigate the idea of a diverging diffusion constant.) Note 
that a slightly different X —>■ oo limit (where the direction of the initial move of the particle 
is maintained until it finds an unoccupied site) was also considered but the product measure 
is still retained. 



In a different context Barma and Ramaswamy |^5[ introduced the 'drop-push' model 



of activated flow involving transport through a series of traps. Each trap can only hold a 
finite number of particles. For the trap depth set equal to one this model is essentially the 



same as the sandpile model of I^Sf discussed above {i.e. it is a zero-range process with some 



infinite rates). In fact the version studied in |7l| is precisely the limit of u{n) — >■ oo for n > 1 



of the totally asymmetric zero-range process described in Section 3.1. A generalisation to 
inhomogeneous traps, and partially asymmetric hopping rates dependent on the occupancy 
of the trap was made in ||7^ and a steady state similar to (11|,1^) demonstrated. 



The zero-range process is also relevant in the context of 1-1-1 dimensional interface growth 
by the step flow mechanism. The interface can be visualised as an ascending staircase of 
terraces. Adatoms land on the terraces and diffuse until they bind to the ascending step. If 
the ratio of deposition rates over diffusion rates tends to zero then the resulting dynamics is 
that a terrace shrinks by one unit (and the adjacent higher terrace grows by one unit) with 
a rate proportional to the size of the terrace. Thus the terrace lengths are equivalent to 
the site occupancies of an asymmetric zero-range process that was discussed in Section 3.1. 
The equivalence of zero range processes to a general class of step flow models is discussed 
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Finally we note that the repton model of gel-electrophoresis [THI studied in the case of 



periodic boundary conditions by |£2| is equivalent to an inhomogeneous zero-range process. 
In this case, the particles of the zero-range process represent the excess stored length of a 
polymer which diffuses along the tube of the polymer. The sites in the zero-range process 
represent the segments of the polymer tube and the inhomogeneities in site hopping rates 
reflect the shape of the polymer tube. 



7 Conclusion 

In this work the aims were to give an overview of the area of phase transitions and ordering 
in one-dimensional systems and also to analyse in some detail a particularly simple model, 
the zero-range process. In section 2 several features were identifled which could lead to the 
anomalous behaviour of ordering and phase transitions in equilibrium systems: long-range 
interactions; zero temperature; unbounded local variable. For nonequilibrium systems some 
concepts which may be important emerged: conserved order parameter; drive; forbidden 
microscopic transitions. 

The simplicity of the zero-range process allowed us to analyse the steady state of the 
model in detail. First we derived the steady state for a general class of zero-range processes 
in Sections 3 and 4. We then analysed the condensation transitions that can occur. On an 
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inhomogeneous system the condensation is very reminiscent of Bose-Einstein condensation. 
For it to occur requires certain conditions to hold on the distribution of hopping rates. In the 
homogeneous system the condensation corresponds to a spontaneous symmetry breaking, 
since an arbitrary site is selected to hold the condensate. The condition for it to occur is 
that the hopping rate dependence on the site occupancy decays sufficiently slowly. It was 
also shown that when the condition for condensation does not hold, one can still observe 
very sharp crossover behaviour and apparent condensation on a finite system 

An interesting possibility that was explored was that of the existence of an effective 
energy function. We saw that any steady state of the form ( pT| , |T7D can be obtained from 
a process obeying detailed balance. However when the effective energy is inferred for cases 
where phase transition occurs (as was carried out for an explicit example in section 5.3) 
we find that it contains long-range interactions. Thus the condensation transition can be 
rationalised within the equilibrium framework. 

Moreover in the zero-range process the existence of a drive or preferred direction, pro- 
ducing a conserved particle current, is not essential for the occurrence of a condensation 
transition. What does appear necessary, however, is the conservation of particles. The fixed 
number of particles implies implies the introduction of a fugacity z through (^2|), which in 
turn controls the condensation transition. As we saw the fugacity gives the hopping rate 
out of a site (referred to as speed in section 5) which is a conserved quantity. 

On the other hand, for other models the presence of a preferred direction and conserved 
current does seem crucial for the existence of phase transitions. For example, the asymmetric 



exclusion process defined in Section 2.2 has non-trivial phase behaviour but the undriven 



version (symmetric exclusion) does not. 

In summary, although a general theoretical framework for the description of phase tran- 
sitions in one dimensional systems is not yet available, we hope that the issues and models 
discussed in the present paper serve to show that our understanding is developing. 
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